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The influence of the atomic-scale inhomogeneities of the pairing interaction strength on the su- 
perconducting order parameter and the conductance spectra measurable by STM is studied in the 
framework of weak-coupling BCS-like theory for two-dimensional lattice model. First of all, it is 
found that the inhomogeneity having the form of atomic-scale regions of enhanced pair interaction 
increases the ratio of the local low-temperature gap in differential conductance spectra to the local 
temperature of vanishing the gap 2Ag/Tp. Even in the framework of mean-field treatment this ratio 
is shown to be larger than the one corresponding to the homogeneous case. It is shown that the 
effect of thermal phase fluctuations of the superconducting order parameter can further increase 
this ratio. Taking them into account in the framework of a toy model we obtained the ratio 2Ag/Tp 
to be ~ 7 — 8. It is found that the additional atomic-scale hopping element disorder and weak po- 
tential scatterers, which can also take place in cuprate materials, have no considerable effect on the 
statistical properties of the system, including the distribution of the gaps, Tp and the ratio 2Ag/Tp. 
The second consequence of the atomic-scale order parameter inhomogeneity is the anticorrelation 
between the low-temperature gap and the high-temperature zero-bias conductance. The obtained 
results could bear a relation to recent STM measurements. 

PACS numbers; 74.45.-|-c, 74.50.-|-r 



I. INTRODUCTION 

Nanoscale inhomogeneities have been widely ob- 
served in the high-temperature superconductor 
Bi2Sr2CaCu208+x (BSCCO) and have generated 
intense interesti'^'^'^'^'^. In particular, the spectral gap 
in the local density of states (LDOS) has been investi- 
gated by scanning tunneling microscopy (STM). It was 
found that deeply in the superconducting state the low- 
temperature gap varies by a factor of 2 over distances 

o 

of 20 — 30 A. Recently the existence of very similar 
picture of inhomogeneities in La^-xStxCuO/^ has been 
also reported^. Several scenarios have been proposed to 
understanding this electronic inhomogeneity. First of 
all, it was speculated that poorly screened electrostatic 
potentials of the dopant atoms vary a doping concentra- 
tion locally, giving rise to the gap modulat ions^i^ i i o , 1 1 ^ 
Alternatively, these inhomogeneities are associated with 
a competing order parameter, such as spini^ii^iii^ and 
orbital antiferromagnetism^^^, charge density wave^^ or 
pair density wave^^. Further, the positive correlations 
between the inhomogeneities and positions of the dopant 
atoms have been observed by STM on the optimally 
doped BSCCO^. After that it was proposed by Nunner 
et al. in Ref. ^ that the dopant atoms modulate the 
pairing interaction locally on the atomic scale. The 
LDOS calculated in the framework of this model is 
in good agreement with the key characteristics of the 
experimental spectra. However, there is an alternative 
picture of the inhomogeneity origin based on the local 
variation of the doping concentration, which explains 
why the correlations are rather weakli. 

On the other hand, in the high- Tc superconductors a 
partial gap in the LDOS exists for a range of tempera- 



tures above T^. There is no consensus up to now if this 
gap is due to pairing without phase coherence, a compet- 
ing order or proximity to the Mott state2£i2ii^i2^. The in- 
homogeneities described above complicate the situation. 
Only very recently the spatially resolved STM measure- 
ments of gap formation in BSCCO samples with different 
Tc corresponding to hole concentrations from 0.12 to 0.22 
were performed^^. For a range of doping from 0.16 to 
0.22 they have found that gaps nucleate in nanoscale re- 
gions above Tc and proliferate as the temperature is low- 
ered, evolving to the spatial distribution of gap values in 
the superconducting state. It was observed experimen- 
tally that overdoped and optimally doped samples have 
identical gap-temperature scaling ratios. Taking into ac- 
count this finding together with the fact that in the over- 
doped samples pseudogap effects are believed to be weak 
or absent and consistency of the low-temperature spec- 
tra with a d-wave superconducting gap. Gomes et.al^ 
have interpreted the gaps above Tc as those associated 
with pairing. Despite the inhomogeneity, every pairing 
gap develops locally at the temperature Tp following the 
relation 2Ag/Tp = 7.9 ±0.5 in wide range of doping from 
overdoped to optimally doped samples. This local pair- 
ing criterion seems to fail only in underdoped samples 
and DOS indicates the presence of another phenomenon, 
possibly unrelated to pairing. 

It is well-known that in the framework of weak- 
coupling BCS theory with homogeneous pairing ampli- 
tude the ratio 2Ag/Tc is 3.5 for s-wave superconductors 
and approximately is in the interval 4.3 — 4.5 for d-wave 
superconductors (this value depends slightly on the par- 
ticular tight-binding parameters). The ratio 2Ag/Tp ^ 
4.7 — 5.2 for d-wave case is a bit higher due to thermal 
smearing of the measured dl/dV curves. It is worth to 
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note that in the framework of strong-coupling theory the 
ratio l^gjTc is in the range 3.5 — 5 for s-wave supercon- 
ductors, but becomes dependent on Ag and can reach the 
values ~ 10 for d-wave pairing case^^. 

In the present paper we show that if the superconduct- 
ing order parameter (OP) is modulated locally on the 
atomic scale, the ratio 2Ag/Tp strongly increases in the 
framework of conventional weak-coupling theory. In ad- 
dition, this inhomogeneity inevitably leads to the anticor- 
relation between the low-temperature gap and the high- 
temperature zero-bias conductance. This findings resem- 
ble the results observed in recent STM experiments2ii2^. 
The most natural way to obtain the atomic-scale modula- 
tions of the OP is to assume the atomic-scale modulations 
of the pairing interaction strength, as it was proposed by 
Nunner et. a&^. We do not mean any particular mech- 
anism of pairing. Our analysis is phcnomcnological and 
the conclusions are independent on the underlying pair- 
ing mechanism and the particular course of the local pair 
interaction modulations. Here we only focus on the ef- 
fect of the inhomogeneity on some observable properties 
of cuprate superconductors. As it is discussed below, if 
this inhomogeneity has the form of atomic-scale regions 
of enhanced pair interaction, the ratio 2Ag/Tp increases 
even in the framework of mean-field treatment. However, 
the effect only takes place if the total area of these re- 
gions is less than the area occupied by the background 
interaction. The last condition is essential. If this is not 
the case, that is the total area of the enhanced pairing 
regions is of the order of or larger than the area of the 
background pairing, the discussed ratio is reduced or, at 
least, remains equal to the homogeneous one. It was 
demonstrated in a number of papers. In particular, the 
attractive Hubbard model with inhomogeneous pairing 
amplitude was studied^ and it was shown that the su- 
perconducting critical temperature can be significantly 
increased in such a system as compared to a uniform 
system corresponding to the pairing interaction averaged 
over the system. The zero temperature superconduct- 
ing order parameter increases also, but in a less degree 
than the critical temperature, what results in reducing 
the ratio 2Ag/Tc. The reduction of the ratio of the en- 
ergy gap to the critical temperature due to the inhomo- 
geneity of coupling constant was also obtained for dirty 
s-wave superconductors^. The potential disorder, as was 
demonstrated-^^, also diminishes this quantity because 
the superconducting order parameter is suppressed on 
the distance of the order of the coherence length around 
an impurity and, consequently, the total area of the en- 
hanced order parameter regions dominates. 

Further, it is physically reasonable that the phase of 
the superconducting order parameter should fluctuate 
from one region of enhanced pairing amplitude to another 
in such an inhomogeneous situation, especially for short 
coherence length cuprate superconductors. In particu- 
lar, the state with antiferromagnetic and nano-scale su- 
perconducting domains exhibiting randomly distributed 
phases was recently proposed^ to account for the for- 



mation of the Fermi arcs, observed in the pseudogap 
phase of the underdoped cuprate a^^'^^'^^ . The existence 
of the thermal phase fluctuations of the superconducting 
clustered state in disordered s-wave superconductors and 
their role in the superconductor-insulator transition was 
demonstrated"^^. In the present paper we show that ther- 
mal phase fluctuations can significantly suppress the tem- 
perature Tp and, consequently, increase the ratio 2Ag/Tp. 
Even in the framework of very simple model we study 
here this quantity reaches the value '^7 — 8, comparable 
to the experimentally observed^l. 

The influence of the additional disorder such as the 
atomic-scale inhomogeneities of the hopping matrix ele- 
ments and weak potential scatterers is also considered. 
It is shown that while they can affect the shape of low- 
temperature LDOS in the system, the properties of inter- 
est: gap's distribution, Tp's distribution and the distribu- 
tion of the ratio 2Ag/Tp remain qualitatively unchanged. 

The paper is organized as follows. Sec. [Til is devoted 
to the detailed mean-field consideration of the problem. 
The model mean-field Hamiltonian we use and the out- 
line of the T-matrix method are introduced in Sec. Ill Al 
The single pair interaction perturbation is studied and 
the physical reasons for the enhancement of the ratio 
2Ag/Tp are discussed in Sec. Ill Bl Sec. Ill CI is devoted 
to the mean-field treatment of the interaction between 
many OP scatterers. The effect of additional weak po- 
tential and hopping clement inhomogeneities is investi- 
gated in Sec. IIIDI The influence of the thermal phase 
fluctuations on the quantities under consideration is dis- 
cussed in Sec. IIIII In Sec. IIVI it is demonstrated that 
the atomic-scale OP inhomogeneity leads to the anticor- 
relation between the local low-temperature gap and the 
value of zero-bias conductance at higher temperatures. 
The conclusions are presented in Sec. |Vl 

II. MEAN-FIELD TREATMENT 

A. Model and method 

We consider the following Hamiltonian on a square lat- 
tice 

H = - Ujcl^Cj^-Y^ Hcl^da+Y^ (Ay4^c]j^ + h.c^ 

ij,a [ij) 

(1) 

where Cia{c\^) stands for an electron annihilation (cre- 
ation) operator at site i with spin tr. 'Yliij indicates sum- 
mation over neighboring sites, while denotes the 
summation over nearest neighbors, ty is the hopping in- 
tegral between sites i and j. We set to be i = 1 for the 
nearest-neighbor hopping and all the energies are mea- 
sured in units of t throughout the paper. The nearest- 
neighbor d-wave order parameter should be determined 
self-consistently: Aij = —gij{ciiCj^ — Cj|Ci|). 

In order to analyze the inhomogeneous pairing correla- 
tions in the framework of the mean-field Hamiltonian ([T]) 
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we exploit the fully self-consistent T-matrix technique for 
Gor'kov Green's functions. Starting from the Gor'kov 
equations the normal and anomalous Green's functions 
are expressed in terms of the homogeneous background 
Green's functions G?, and the T-matrix, which contains 
all the inhomogeneities. Then the full Green's function, 
which depends on two space indices i and j, takes the 
form 

Gio — G^ij + ^ G^k'i'kmGm] ■ (2) 

Here fkm = ~J2i^'^~^)knVnm, Mkm = Skm + 

n 

J^^kn^nm- Vkm IS the perturbation matrix including 

n 

all the inhomogeneities. All Green's functions and T- 
matrices are 4x4 matrices in the direct product of spin 
and particle-hole spaces, what indicated by the symbol 

fi and Ui are Pauli matrices in particle-hole and spin 
spaces, respectively. The summation is taken over all the 
sites, where the OP = A°„j -f (5Afe,„ differs from the 
background value A°^ . A°^ is assumed to be of d-wave 
type, that is A? i ^ = — A° - = A°. a and b are basis 
vectors of the square lattice. We set the lattice constant 
a to be equal to unity. 

OP is to be calculated from the self-consistency equa- 
tion 

Ay = 5yr^Tr4 [f_io-yGij(e„)] : (3) 

where = (fx — ify)/2 and e„ is the fermionic Mat- 
subara frequency. The energy cutoff |£„| < 4 is used in 
Eq. ^ for the technical convenience. We have checked 
that the particular value of the cutoff energy (if it is 
considerably larger than the maximal value of the OP) 
does not change qualitatively the results and only slightly 
varies the OP quantitative values. Eqs. ^ and ^ al- 
low us to find the OP A^j numerically. Making use of 
the outlined technique we consider a square of the size 
n X n, immersed into an infinite lattice carrying the ho- 
mogeneous OP A? , as an inhomogeneity described by 
T-matrix. 



B. Single perturbation 

In general, the T-matrix can describe inhomogeneities 
of all the parameters entering the equations: hopping 
elements, chemical potential or self-energies. For the 
moment we only focus on the off-diagonal self-energy 
inhomogeneity and therefore Vkm = SAkmio'2if2. We 
begin by considering a single perturbation of the pair- 
ing interaction. Two different models for the individ- 
ual pairing perturbation are represented in the paper. 
Model (a) describes the region of the enhanced pair in- 
teraction emanating from a site i by Yukawa-type po- 
tential gij = gb + 5gfij, where gb corresponds to a 
background interaction and fij = J2s^^Pi~''^i3s / ^) /''^ijs- 



Tijs = Y + is the distance between the cen- 
ter of the bond connecting the sites i and j and the 
source of the pairing interaction perturbation beyond 
the plane. Here z is the distance between the Cu02 
plane and the source of the perturbation and r^^l^^'J de- 
notes the distance in the plane. Model (b) supposes that 
the pair interaction is considerably enhanced on a pla- 
quette and has a long weak tail beyond this plaquette: 
Qij = gb + ui if the bond ij belongs to the chosen plaque- 
tte and gij = gb + ^2/(^0 ^ +2:^) if the bond ij is beyond 
this plaquette. Here, ro.y is the distance between the 
center of the chosen plaquette and the center of the bond 
ij. By considering a number of different models for the 
single perturbation we have checked that the results do 
not depend qualitatively on the particular shape of the 
perturbation and are only controlled by its effective width 
and height. The tight binding model parameters taken 
to obtain the normal quasiparticle dispersion are the fol- 
lowing: the next-nearest neighbor hopping t' = —0.3 and 
the chemical potential fi = —1 correspond to the nor- 
mal state background of model (a) and t' = —0.35 and 
fi — —0.8 describe the normal background of model (b). 
Both sets of the parameters give approximately the same 
low-energy dispersion and adjust to qualitatively repro- 
duce experimentally measured Fermi surface of BSCCO 
near optimal doping. However, they result in quite dif- 
ferent quasiparticle dispersions for the energies, which 
have the absolute value of the order of superconducting 
gap, and, in particular, differ by the energy location of 
the normal state van Hove singularity SyH = — M + 4^'. 
This fact leads to essential difference in the shape of the 
LDOS curves, but does not influence qualitatively the 
statistical properties considered here: gap (A^) distri- 
bution, the local temperature of vanishing the gap (Tp) 
distribution and the ratio 2Ag/Tp. We comment on the 
role of the particular set of tight binding model param- 
eters in more detail below. Making use of the T-matrix 
approach outlined above we calculate the superconduct- 
ing order parameter, local density of states and Tp for the 
two inhomogeneous patterns with a single perturbation 
described by models (a) and (b), respectively. 

The resulting space distributions of low-temperature 
site-averaged OP A^ = (A^^i+a + Aj,i_a + |A. .^^| + 
|A^ j_;;|)/4 are presented in Figs.fTJa) andl^^a) for models 
(a) and (b), respectively. 

STM technique measures the local differential conduc- 
tance, and the thermally smeared LDOS, described by 
the expression 

00 

dl/dV^ J ds^i^^^^p,{e,T) , (4) 

— CSO 

can be extracted from these measurements. We set |e| = 
1 throughout the paper. In the above expression 

p,(£,r) = -(l/^)ImG,^(£,r) (5) 

is the local density of states, /(e) is Fermi distribution 
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FIG. 1: (a) Low-temperature (T = 0.03) distribution of the 
site-averaged OP. (b) Low-temperature (T = 0.03) gap map. 
(c) Tp— map. (d) The space profiles of the ratio 2Ag/Tp (solid 
line) and the site-averaged OP (dashed line) calculated along 
the horizontal line drawn though the center of the perturba- 
tion in panel (a). The left (right) vertical axis corresponds 
to the ratio 2Ag{i)/Tp{i) (superconducting order parameter). 
All the pictures correspond to the single perturbation de- 
scribed by model (a) (see text) with gi, — 0.51, Sg = 2.09, 
A = 1.5 and z — 1.5. The parameters A and z are measured 
in units of the lattice constant a. 



function and V is voltage applied between the STM tip 
and the sample. For the homogeneous situation the 
distance between the coherence peaks 2Ag in the low- 
temperature conductance spectra is directly connected 
to the low-temperature OP by the simple relation Ag = 
aAo(T), where the parameter a is slightly model depen- 
dent. In our case a = 4.0 for model (a) and a = 3.6 for 
model (b). 

For an inhomogeneous system, where the characteristic 
size of the patch < ^, there is no any direct simple rela- 
tionship between the local OP and the local gap A^. One 
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FIG. 2: The panels (a)-(d) represent the same quantities as 
in Fig. [T] but correspond to the single perturbation described 
by model (b) (see text) with gb — 0.59, ui = 1.00, ii2 = 0.25 
and z — 0.5. 



can only conclude from the numerical calculations^ii^ 
that for regions, where the OP is enhanced from the 
background, the gap gets wider and the peak height is 
suppressed compared to the average value. It is worth 
noting that, unlike the homogeneous situation, this peak 
does not always represent the maximal superconducting 
gap on the Fermi surface, but in some cases originates 
from the spectral weight transfer from the nearby van 
Hove singularity due to the Andreev scattering processes. 
Otherwise, if the OP in a cluster is less than that one in 
the background, the narrow and high Andreev resonant 
peaks develop in the cluster region resulting in dimin- 
ishing of the gap region. For this reason we investigate 
not only the OP distribution, but also the experimen- 
tally measurable thermally smeared LDOS, which has a 
maximum at V = Ag . For models (a) and (b) of the sin- 
gle perturbation the corresponding low-temperature gap 
maps are represented in Figs.fjb) andlD^b), respectively. 
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Experimentally^'* the temperature Tp{i) of the gap dis- 
appearing for the particular location in the sample has 
been determined using the criterion dI/dV{V = 0) > 
dl/dV (for all > 0). Using the above criterion we cal- 
culated the distribution Tp{i) from the thermally smeared 
LDOS curves. The corresponding Tp-maps for the single 
perturbations described by models (a) and (b) are rep- 
resented in Fig. [IJc) and[2Ic), respectively. It can be 
seen from Figs. [1] and [2] that while the space distribu- 
tion of the superconducting order parameter mainly fol- 
lows that one of the pairing interaction strength due to 
the local multiplication of the anomalous Green's func- 
tion in the self-consistency equation ([3]) by the coupling 
constant, LDOS and Tp are essentially nonlocal (on the 
atomic scale) quantities with the characteristic size of 
order of determined by the scale, where the Green's 
functions change considerably. This is valid for the both 
considered models, however the non-locality is much less 
pronounced for model (b). The reason is that the pair 
interaction perturbation is more spikier in this case as 
compared to model (a). Due to its very small effective 
width the perturbation (b) cannot contribute consider- 
ably to the Green's function and, consequently, to the 
LDOS anywhere in space except for its own plaquette. 
The reason for enhancement of the gap at this plaque- 
tte is spectral weight transfer from the nearby Van Hove 
singularity, which is discussed in detail later. 

The ratios 2Ag{i)/Tp{i) calculated along the horizon- 
tal lines drawn through the centers of the perturbations 
in Figs, [ija) and ^a) for the model samples considered 
above are plotted in Figs.[l](d) and[2l(d) together with the 
spacial profiles of the superconducting order parameter 
along the same cuts. It is seen that this ratio is con- 
siderably enhanced in comparison with the homogeneous 
situation. For model (a) the area of the enhancement is 
larger that the characteristic size of the superconducting 
order parameter (and coupling constant) perturbation. 
This is due to the above mentioned gap and Tp nonlocal- 
ity and leads to the increasing of the ratio over the entire 
sample in the situation with many off-diagonal scatterers 
in spite of the fact that only smaller part of the sample is 
occupied by the enhanced coupling constant regions. On 
the contrary, for model (b) the two profiles practically co- 
incide. Note that the maximal ratio is very high 14) 
for this model. The physical reasons for this fact will 
be given later. The averaged ratio {2Ag{i)/Tp{i)) cal- 
culated over the region of the size, which approximately 
corresponds to the average distance between the centers 
of enhanced pairing regions in the situation with many 
individual perturbations considered below, is equal to 7.6 
and 8.5 for models (a) and (b), respectively. 

There are two main physical reasons for the enhance- 
ment of the ratio 2Ag/Tp for the single perturbation con- 
sidered above. First of all, small as compared to super- 
conducting coherence length perturbations cannot main- 
tain superconductivity by themselves and only do this 
due to the superconductivity in the bulk. Although the 
value of the zero-temperature OP strongly increases when 



the height of the perturbation grows, at finite tempera- 
tures where the bulk OP vanishes, the pairing correla- 
tions in the small area go to zero abruptly. To illustrate 
this the dependence of the site-averaged OP on tempera- 
ture is presented in FiglS] Panels (a) and (b) correspond 
to models (a) and (b) described above. The bottom thin 
line represents the temperature dependence of the back- 
ground order parameter generated by the coupling con- 
stant gb- The bold solid curves correspond to the tem- 
perature dependence of the site-averaged order parame- 
ter for the sites marked by the appropriate numbers in 
Figs. [TJa) andOJa), respectively. Although the super- 
conducting order parameter vanishes at the critical tem- 
perature denoted by over the entire inhomogeneous 
sample simultaneously, as it should be in the mean-field 
approximation, the pairing correlations in the small area 
reduce abruptly near the temperature of the vanishing of 
the bulk OP corresponding to the coupling constant gb- 
Then only the non BCS-like tails extend up to T^. In 
other words, the OP inside the area of the enhanced pair 
interaction disappear with temperature more faster than 
the bulk OP having the same low-temperature value. To 
demonstrate this the temperature dependence of the bulk 
OP having the same low-temperature value as in the cen- 
ter of the perturbation (curves marked by 1) is depicted 
by the upper thin lines in Figs, ^a.) and (b). Note that 
for model (b) the abrupt suppression of the OP near the 
background critical temperature Tb is seen more clearly 
due to smaller effective size of the perturbation. 
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FIG. 3: Dependence of the site-averaged OP on temperature 
for the single perturbation. The thin bottom line corresponds 
to bulk order parameter generated by the background cou- 
phng constant. The bold solid curves represent the tempera- 
ture dependence of the superconducting order parameter for 
the sites marked by the appropriate numbers in Fig. [TJa) and 
[^a). The upper thin line is the temperature dependence of 
the bulk OP with the same zero-temperature value as in the 
center of the perturbation (curve marked by 1). The panels 
(a) and (b) correspond to the models (a) and (b), respectively. 

The temperatures Tp of vanishing the gap for the given 
space locations are marked by the filled circles on the cor- 
responding A(r) curves in Figs, ^a) and (b). It is seen 
that even for the homogeneous case Tp is slightly lower 
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than the critical temperature (see the bottom thin line) 
due to the thermal smearing of the measured differen- 
tial conductance according to the formula ([4]). Conse- 
quently, the bulk ratio 2Ag/Tp is a bit higher than the 
ratio 2Ag/Tt,. These ratios are dependent slightly on the 
particular coupling constant and tight-binding parame- 
ters. For model (a) we consider they have the following 
values: (2Ag/Tb)(") = 4.4 and {2Ag/TpY'''> = 5.0, while 
for model (b) {2Ag/Tb)'^''^ = 4.2 and (2Ag/Tp)(^) = 4.7. 
For the inhomogeneous situation with the single pertur- 
bation the difference between and Tp (Tp is now posi- 
tion dependent) is much more pronounced. This results 
from the discussed above non BCS-like behavior of the 
OP on temperature: for the long OP tails, taking place 
for the temperatures higher than the background criti- 
cal temperature, the gap is too small and easily smeared 
out by high enough temperature. Although at the cen- 
ter of the perturbation the value of the superconducting 
OP is still high enough at the temperatures of the order 
of Tp, the area occupied by this enhanced OP is small 
as compared to ^g. Consequently, the amplitude of the 
corresponding peak in the LDOS is low and again can 
be easy washed out by temperature. These reasons work 
more efficiently for more spikier perturbations. This is 
clearly seen in FiglS] where the difference between Tp 
and Tx is more pronounced for model (b). 

When the size of the pairing interaction perturba- 
tion increases, T^ grows considerably, while the zero- 
temperature value of the superconducting OP in the per- 
turbation region rises only slightly. Naturally, when the 
size of the enhanced pairing area becomes of the order of 
a few superconducting coherence lengths, T^ — > Tf, cor- 
responding to the coupling constant gh + 6g and the bulk 
value of 2Ag/Tp should be restored at the center of the 
cluster. For s — wave pairing case and cluster width com- 
parable to or larger than the reduction of the critical 
temperature of the cluster having the pairing potential 
gi (in comparison to the bulk value corresponding to gi) 
due to the proximity to the background with the pair- 
ing potential go < gi has been studied in Ref.'^^ in the 
framework of quasiclassical Usadel equations. As it is 
seen from the above picture, the OP inhomogeneity of 
the size comparable to or larger than cannot give rise 
to large enough ratios 2Ag/Tp and more spikier atomic- 
scale inhomogeneities are more favorable. 

The second reason for the ratio 2Ag/Tp to get larger 
is the following. As it was already mentioned above, in 
the inhomogeneous situation with an atomic-scale sin- 
gle perturbation of the large enough strength the main 
peak in LDOS does not represent the maximal local 
superconducting gap on the Fermi surface, but origi- 
nates from the spectral weight transfer from the nearby 
van Hove singularity due to the Andreev scattering pro- 
cesses. For the tight-binding parameters we consider the 
van Hove singularity energy in the homogeneous sys- 
tem is approximately described by the formula e^H = 
— + where A^ax is an effective max- 

imal superconducting gap on the Fermi surface. The 



processes of Andreev scattering from the inhomogene- 
ity result in appearing of the symmetrical with respect 
to Fermi energy peak in LDOS. The examples of low- 
temperature LDOS curves exhibiting such a feature are 
shown in Fig. IH Panels (a) and (b) correspond to models 
(a) and (b), respectively. It is seen that for model (a), 
where the normal state van Hove singularity is very close 
to the Fermi energy (and, consequently, the van Hove sin- 
gularity peak in the superconducting state is close to the 
superconducting coherence peak), the above mentioned 
transfer of the spectral weight leads to small enough in- 
crease of the gap Ag . The shape of the resulting peak is 
close to the shape of a superconducting coherence peak 
because of partial overlapping of the van Hove singularity 
and superconducting coherence peak in the bulk. At the 
same time for model (b), where the van Hove singular- 
ity is more distinct from the superconducting coherence 
peak in the bulk, the clearly seen spectral weight transfer 
results in considerable increase of the gap Ag. However, 
this mechanism of a gap enhancement is very local, as 
it is seen in Fig. [H^b). This is the reason for appearing 
of very high gap just at the plaquette occupied by the 
perturbation in Fig. [2Kb). However, the shape of the re- 
sulting peak is strongly distorted. It is worth to note that 
the LDOS curves represented here are related to the sin- 
gle perturbation. For the many individual perturbations 
discussed below the curves have more resemblance to the 
experimentally measured LDOS spectra. 
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FIG. 4: Low-temperature (T = 0.03) LDOS curves taken for 
a number of sites from the center of the perturbation (black 
solid curve) to the background (gray line). The curves cor- 
respond to the sites marked by the appropriate numbers in 
Fig. [ija) and Fig. [2ja). Panels (a) and (b) correspond to 
models (a) and (b), respectively. 



C. many OP scatterers 

Now we turn to discussion of more realistic situation, 
when many pair interaction scatterers are present in the 
sample. We have considered a 22 x 22 sites square as a 
perturbation described by the T-matrix. We again repre- 



sent here two models. In the both models individual per- 
turbations are randomly distributed in the square with 
the concentration n — 0.07 (here concentration means 
the number of individual scatterers divided by the to- 
tal number of sites in the square). In model (A) the 
individual perturbations have the same shape as in the 
single perturbation model (a). Model (B) corresponds 
to the individual perturbations described by model (b). 
The tight-binding parameters describing the normal state 
quasiparticle dispersion coincide for models (a) and (A) 
so as for models (b) and (B) , respectively. The numerical 
values of the parameters, characterizing the individual 
scatterers, are taken to give the same low-temperature 
values of superconducting OP in the center of the per- 
turbation and in the background as for the corresponding 
single perturbation model (see captions to Figs. [H] and [S] 
for the particular values). 

Since we consider the finite-size square as a pertur- 
bation described by T-matrix, the finite-size effects will 
affect the results. However the appropriate choice of the 
homogeneous superconducting OP A" outside the square 
helps us to minimize the influence of the square bound- 
ary. By considering a larger 30 x 30 sites square we have 
checked that if the bulk OP outside the square is taken 
to be equal to the average background value of the OP 
in the square, the error in calculating OP, gap and Tp 
becomes less than 5% at the distance ~ 4 sites from the 
boundary and lessens further towards the center of the 
square. For this reason the 14 x 14 central region of the 
square is taken into account in calculating the statistical 
properties of the system. 

The resulting spatial distributions of low-temperature 
site-averaged OP, low-temperature gap-maps and Tp 
maps are demonstrated in Figs. [5] and [6] for models (A) 
and (B), respectively. It is worth to note that the tem- 
perature evolution of OP spatial distribution was already 
discussed in Ref. [H, so we do not focus here on this as- 
pect. The spatial distributions of low-temperature OP 
are represented in Figs.jHKa) andl^Ka). As it was described 
earlier for the single perturbation, the spatial profile of 
the OP inhomogeneity mainly follows that one of the 
coupling constant, that is the characteristic scale of the 
corresponding inhomogeneity is ~ 2 atomic sites. How- 
ever, due to discussed above fact that the deviation of 
the Green's function from the bulk value is spread over 
a wider space region than the superconducting OP, the 
gap and Tp inhomogeneities are more smooth and the 
characteristic size of the patches in the gap and Tp maps 
is of the order of 4 — 5 atomic sites for model (A) and ^ 3 
atomic sites for model (B). For model (A) the size of a 
patch is roughly in accordance with the superconducting 
coherence length ^g, while for model (B) it is somewhat 
smaller due to very spiky character of the perturbation. 

The probability distributions to find the particular 
value for 2Ag{i)/Tp{i) in our model samples are plotted 
in Figs.[5Jd) andlHId). It is seen that the distribution for 
model (A) is quite narrow. The reason is that the spatial 
profile of the ratio 2Ag/Tp is broadened and flattened as 




FIG. 5: (a) Low-temperature (T = 0.05) distribution of the 
site-averaged OP for model (A), (b) Low-temperature gap 
map. (c) Tp— map. (d) Histogram of 2Ag{i)/Tp{i) values. 
The parameters describing a perturbation are the following: 
Qb = 0.05, Sg = 2.72, A = z = 1.5. 



compared to OP spatial profile, as it was demonstrated 
in Fig. [Hd) and discussed in the context of the single 
perturbation. For model (B) the distribution is wider 
and exhibits a long tail of very high ratios for a small 
part of sites. The reason for this fact was discussed in 
the context of the single perturbation. 

The average values of the ratio 2Ag{i)/Tp{i) for models 
(A) and (B) are 5.4 and 5.6, respectively. These values 
are considerably less than the ones for the appropriate 
single perturbations, described by models (a) and (b) (as 
it was discussed earlier, for the single perturbation the ra- 
tio {2Ag{i)/Tp{i)) averaged over the region ~ P, where 
I corresponds to the average distance between the cen- 
ters of enhanced pairing regions for the appropriate many 
scatterers model, equals to 7.6 and 8.5 for models (a) and 
(b), respectively). Although the mean-field ratios for 
many scatterers models exceed the corresponding bulk 
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FIG. 6: (a) Low-temperature (T = 0.03) distribution of 
the site-averaged OP for model (B). (b) Low-temperature 
gap map. (c) Tp— map. (d) Histogram of 2Ag{i)/Tp{i) val- 
ues. The individual perturbation and the background are 
described by the following parameters: gt = 0.25, ui — 1.16, 
U2 = 0.25 and z — 0.5. 



ratios 2A^'''7r/''"' = 5.2 and 2Af ""/Tp^''^" = 4.8, this 
increase is quite small as compared to the single pertur- 
bation, especially for model (A). Here the homogeneous 



ratios 2A 



A{B},av ,rpA{B),av 



are calculated for the coupling 



constants, which give the same zero-temperature OP, as 
the average zero-temperature OP in samples (A) and 
(B). The reason for such a crucial reduction of the ra- 
tio 2/S.g/Tp in the presence of many OP scatterers is 
the proximity effect. The individual perturbations essen- 
tially interact one with another for the considered con- 
centrations when the averaged distance between them is 
less then the superconducting coherence length. In the 
framework of the mean-field treatment the phase of the 
OP is the same over the entire sample. For this reason 
the individual perturbations raise the pairing correlations 
in the background and maintain their own OP even for 



high enough temperatures in comparison with the back- 
ground critical temperature. This is quite different from 
the picture of the single perturbation, where the pairing 
correlations sharply weaken for the temperatures higher 
than Tfe. 
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FIG. 7: The dependence of the site-averaged OP on tem- 
perature for the sites marked by the appropriate numbers in 
Figs. [Ha) and|6l^a) (black solid lines). The gray lines are the 
temperature dependence of the bulk OP having the same zero- 
temperature value as in the center of the chosen perturbation 
(upper gray curve), at the chosen site belonging to the back- 
ground (bottom gray curve) and as the averaged OP (middle 
gray curve). The dotted line represents the temperature de- 
pendence of the OP, averaged over the entire sample. Panels 
(a) and (b) correspond to models (A) and (B), respectively. 

The discussed behavior of the superconducting OP is 
depicted in Fig. [T] Panels (a) and (b) correspond to 
models (A) and (B), respectively. The black solid lines 
represent the dependence of the site-averaged OP on tem- 
perature for the sites marked by the appropriate numbers 
in Figs, in^a) and[6Ua). They correspond to the charac- 
teristic OP behavior along the line from the center of the 
perturbation to the background (from top to bottom). 
The gray lines illustrate the temperature dependence of 
the bulk OP having the same zero-temperature value 
as in the center of the chosen perturbation (upper gray 
curve), at the chosen site belonging to the background 
(bottom gray curve) and as the averaged OP (middle 
gray curve). The dashed line represents the temperature 
dependence of the OP, averaged over the entire sample. 
Although the curves still exhibit non-BCS shape, there is 
no sharp suppression of the superconducting correlations 
in the vicinity of the background critical temperature un- 
like the single perturbation case. Now, roughly speak- 
ing, it is the critical temperature of the bulk OP, hav- 
ing the same low-temperature value as the averaged OP 
(middle gray curve), that plays part of the background 
critical temperature for the single perturbation: if the 
OP at a given location is less than the averaged value, 
the superconducting correlations for a finite temperature 
are higher than the bulk ones corresponding to the same 
zero-temperature OP, and only for the locations, where 
the OP exceeds the average value, the finite-temperature 
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superconducting correlations are suppressed as compared 
to the bulk behavior. This leads to the fact that the val- 
ues of Tp{i) are higher than for the single perturbation 
case, what, in turn, results in considerable decrease of the 
ratio 2Ag/Tp. It is worth noting that in our mean- field 
treatment we neglect OP phase. However, it is physically 
reasonable that the phase of the superconducting order 
parameter should fluctuate from one region of enhanced 
pairing amplitude to another in such an inhomogeneous 
situation, as it was already mentioned in the introduc- 
tion. We discuss the influence of the thermal phase fluc- 
tuations on the above results in the next section. 

As it was mentioned above, the difference between the 
averaged ratio {2Ag{i)/Tp{i)) for model (B) and the cor- 
responding homogeneous value is higher than between 
the same quantities in model (A) . In addition, the proba- 
bility distribution to find the particular value of the ratio 
has a long tail of very high values for this model. The ra- 
tio is greatly enhanced for a small part of the sites due to 
the fact that the mechanism of the spectral weight trans- 
fer from the van Hove singularity makes the gap consid- 
erably larger for this model, where the van Hove singular- 
ity and the bulk coherence peak are quite distinct, but it 
only works in the close vicinity of an individual perturba- 
tion. To illustrate this the LDOS and thermally smeared 
LDOS calculated according to the expression ^ at some 
typical sites are shown in Figs.[51[a),(b) for model (A) and 
Figs. [8Uc),(d) for model (B), respectively. The left col- 
umn represents the bare LDOS, while the right one cor- 
responds to the thermally-smeared quantity. Please note 
that the temperatures we use to calculate the thermally 
smeared LDOS are deeply in the superconducting state, 
where the superconducting OP practically does not differ 
from its zcro-tcmpcrature value. 

It is seen from the figures that for the most part of 
sites the bare LDOS exhibits two peaks (and also a num- 
ber of kinks). The peak at smaller energies represents a 
superconducting coherence peak corresponding to an ef- 
fective OP. This effective OP neither coincides with the 
background value, as it was for the single perturbation, 
nor represents the local OP value. Apparently, its value 
is close to the average over some area for the most part of 
the samples, however this will be discussed in more detail 
elsewhere^. The second peak at higher energies is usu- 
ally weaker (except for the small part of sites in the sam- 
ple (B)) and associated with the spectral weight transfer 
from the van Hove singularity due to Andreev scattering. 
For the both models the spectral weight transfer is clearly 
seen in the bare LDOS, however the two peaks overlap 
for model (A) due to the fact that the normal state van 
Hove singularity is quite close to the Fermi surface in this 
case. As a result the effective superconducting coherence 
peak and the one transferred from the van Hove singu- 
larity merge in the thermally smeared LDOS even for a 
low temperature, as it is seen in Fig[8fb). The resulting 
peak resembles a superconducting coherence peak. 

For model (B) the effective superconducting coherence 
peak and the peak transferred from van Hove singular- 
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FIG. 8: LDOS (left column) and thermally smeared LDOS 
(right column) . For all the figures the black solid curve repre- 
sents LDOS at a given site typical for the background, dashed 
line corresponds to the center of a perturbation and LDOS 
somewhere between them is plotted by the dotted line. Upper 
row represents the results for model (A), thermally smeared 
LDOS is calculated at T = 0.05. Middle row demonstrates 
the results for model (B), while bottom row corresponds to 
model (B'). Thermally smeared LDOS for models (B) and 
(B') is obtained at T = 0.03. 



ity are more distinct. In this case the effective coherence 
peak is higher than the transferred one for a number 
of sites. However, they are of the same order for the 
most part of sites. As a result the gap position is lo- 
cated somewhere between them in the low-temperature 
thermally smeared LDOS. The corresponding curves do 
not exhibit a pronounced peak and resemble the dl/dV 
behavior observed in underdoped samples in contrast to 
model (A) LDOS, which is more similar to optimally and 
overdoped samples. Finally, the transferred peak wins for 
a small number of sites, typically at centers of the indi- 
vidual perturbations. This leads to the sharp increase of 
the gap for these spatial locations and, consequently, to 
the existence of the long tails in the 2Ag/Tp distribution. 



D. effect of weak potential and hopping element 
inhomogeneities 

It is reasonable to assume that the possible causes 
of the pairing interaction inhomogeneities (for example, 
dopant atoms) most probably also give rise to a disor- 
der of the normal state parameters. Therefore, let us 
turn to the discussion of the additional effect of inho- 
mogeneities of tight binding parameters such as chemical 
potential and hopping matrix elements on the measured 
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by STM properties. Here we only focus on a weak in- 
homogeneity of the chemical potential \Sfi\ < t. Strong 
potential scatterers characterizing by \Sfi\ ^ t are well- 
known to give rise to quasiparticle resonant bound states 
in the close vicinity of the potential impurity^^. In real 
BSCCO system only very small part of in plane native 
defects exhibits near zero bias resonances, so they should 
not influence considerably the statistical properties stud- 
ied here, such as the 2/S.g/Tp distribution. 

To illustrate the influence of this additional disorder we 
have chosen our model (B), where the nearest neighbor 
hopping element is taken to be t -I- 6t, the next-nearest 
neighbor equals to t' + 5t' and the chemical potential is 
^ + Sfj, for the plaquettes corresponding to the enhanced 
pairing interaction. Here we take 6t = 0.50, 6t' = 0.15 
and — —0.25. This model is referred to as (B'). 

The LDOS and low-temperature thermally smeared 
LDOS for model (B') are shown in Figs. [5]^e) and (f), 
respectively. As it is seen in the figures the low-energy 
part of the curves is practically not affected by the po- 
tential and hopping matrix element disorder. Only the 
high energy part and, especially the region of the van 
Hove singularity is distorted by this type of scatterers. 
The van Hove singularity become less pronounced under 
the influence of the additional disorder: more widened 
and reduced in height. This is natural because in the 
homogeneous situation the change of tight binding pa- 
rameters strongly shifts the energy location of the van 
Hove singularity. It is worth to note that for the mod- 
els, where the normal state van Hove singularity is more 
close to the Fermi energy (similar to model (A) we con- 
sider), the effect of the diagonal scatterers on the LDOS 
is even weaker, because the van Hove singularity is ab- 
sorbed by the superconducting coherence peak and the 
resulting peak is not qualitatively sensitive to the dis- 
cussed types of disorder, as it was shown in Ref. [H for 
the case of weak potential disorder. 

To investigate the influence of the discussed additional 
disorder on the properties of interest here we compare 
the OP, gap and Tp spatial distributions for models (B) 
and (B'). The correlation between the OP spatial dis- 
tributions in models (B) and (B') is demonstrated in 
Fig. [9{a). The horizontal axis represents the value of 
the superconducting OP in model (B), while the vertical 
axis corresponds to this value for model (B'). Each point 
with the coordinates (A^-^^, A'^ represents the values 
of the superconducting OP for a given site i in models 
(B) and (B'). The line is the linear fit A^-^') = I.OIA^-^) 
to these points. It is seen that the points are described 
by this fit very well, so the additional disorder practically 
docs not affect the superconducting OP. 

At the same time the gap and Tp distributions are in- 
fluenced by the potential and hopping scatterers, as it is 
seen in Figs. [Sib) and (c), respectively. The reason is that 
the gap is (at least partially) determined by the LDOS 
peak originated from the spectral weight transfer from 
the van Hove singularity by Andreev scattering processes, 
while, as it was discussed above, the van Hove singular- 
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FIG. 9: (a) The correlation between the OP spatial distribu- 
tions in models (B) and (B'). Each point with the coordinates 
( A'^' , A'^ ' ) represents the values of the superconducting OP 
for a given site i in models (B) and (B'). The line is the linear 
fit A'^') = I.OIA'-^' to these points, (b) Analo sous correla- 
tion between the low-temperature (T = 0.03) gaps in the con- 
sidered models. The points are fitted by the liner dependence 
A^-^'^ = 0.69A^-^' -I- 0.18. (c) The correlation between Tp. 

The linear fit is T^^'' = 1.42r^^' - 0.06. (d) The histograms 
for 2/S.g/Tp distributions. The upper histogram corresponds 
to model (B'), while the bottom is related to model (B). The 
offset is for clarity. 



ity is affected quite strongly by the potential and hopping 
disorder. Nevertheless, the correlations between the low- 
temperature gaps and Tp in models (B) and (B') can be 
still fitted by a linear dependence, what indicates the fact 
that the corresponding gap maps do not change qualita- 
tively under the influence of such a disorder, but only 
distorted to a certain extent. The average values of the 
low-temperature gap and Tp also differ very slightly for 

models (B) and (B'): (A^^'^) = 0.50 and {T^^'^) = 0.17, 
while (A^^') = 0.47 and (T^^') = 0.16. 

The distributions of the ratio 2Ag/Tp for models (B) 
and (B') are compared in Fig.[5Jd). The average value of 
the ratio is again affected by the potential and hopping 
disorder only slightly: {2Ag/Tp) = 5.8 for model (B') 
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and (2Ag/Tp) = 5.6 for model (B), although the disorder 
results in widening of the distribution. 

At the end of the discussion of the potential and hop- 
ping disorder effect we would like to note that the con- 
sidered here set of parameters is just a representative 
example. We have studied a number of other weak po- 
tential and hopping disorder configurations and found 
that their effect qualitatively the same. Although the 
particular values for the potential and hopping disorder 
parameters and their relationship to the pair disorder 
strength can be only established in the framework of cer- 
tain microscopic models, our analysis indicates that they 
do not qualitatively influence the physics discussed in the 
present paper. 



III. THERMAL PHASE FLUCTUATIONS: A 
TOY MODEL 

The considered above mean-field approximation seems 
to be physically not quite appropriate for studying the 
state, which is inhomogeneous on the atomic scale, espe- 
cially in view of short coherence length in cuprate ma- 
terials. It is reasonable to assume that the phase of the 
superconducting order parameter should fluctuate from 
one region of enhanced pairing amplitude to another 
in such an inhomogeneous situation, what can signifi- 
cantly suppress the temperature Tp and, consequently, 
increase the ratio 2A/Tp. The regular consideration of 
the state with inhomogeneous pairing interaction beyond 
the framework of mean-field approximation, which takes 
into account thermal phase fluctuations, is a separate 
problem. So we postpone it for a future publication. To 
show that the thermal phase fluctuations are indeed sys- 
tematically suppress Tp, in the present paper we study 
their effect in the framework of a toy model. 

While in the mean-field approximation the phase of the 
superconducting OP is the same over the entire sample, 
we assume that it only remains constant over the region 
^ P around an individual perturbation, where I is an av- 
erage distance between the perturbations. The value of 
the phase in the vicinity of a perturbation centered at a 
site i for a given temperature T is set by hand according 
to the formula 7T{0i}, where is a random number 
belonging to the interval [— 7r/2,7r/2] and 7 is a coeffi- 
cient accounting for the strength of the fluctuations. As 
it is seen from the above expression the phases at differ- 
ent perturbation regions are partially correlated for low 
enough temperatures and the variation rises with tem- 
perature modelling the effect of the thermal fluctuations. 
The particular law of increasing the fluctuations with ris- 
ing temperature (linear in the considered model) is not 
very important. Anyway, the main effect of the thermal 
phase fluctuations is to suppress Tp considerably. The 
reason for this is that they partially destroy the prox- 
imity effect between different perturbation regions thus 
shifting the physical properties in the vicinity of a given 
pair scatterer to the limit of an independent single per- 



turbation. 

To demonstrate the discussed above effect in the frame- 
work of our toy model, in Figs. fTUf a) and fTTT a) we plotted 
the part of the sample gapped at a given temperature as 
a function of temperature for models (A) and (B), respec- 
tively. The black solid line in Fig.fTUTa') [HJa)] represents 
the results for the pure mean-field model (A) [(B)] with- 
out phase fluctuations, the dashed lines demonstrate the 
effect of weak enough fluctuations corresponding to 7 = 5 
and the dotted lines are related to more stronger fluctu- 
ations with 7 = 10. It is seen that the increase in the 
fluctuation strength 7 monotonously suppresses Tp. 
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FIG. 10: (a) The part of the sample gapped at a given tem- 
perature as a function of temperature. The black solid line 
corresponds to the mean-field results, the dashed line demon- 
strates the effect of weak enough fluctuations described by 
7 = 5 and the dotted line is related to more stronger fluctu- 
ations with 7 = 10. (b) The influence of the fluctuations on 
the 2Ag{i) /Tp{i) distribution. The bottom histogram is cal- 
culated for the mean-field sample, while the middle one is for 
the weaker fluctuations with 7 = 5 and the upper distribution 
corresponds to 7 = 10. All the results are related to model 
(A). 
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FIG. 11: The same results as in Fig. 1101 but calculated for 
model (B). 

Figs. [TUTb') and [TlT b) represent the influence of the 



12 



7 





5 


10 


(A) 


5.4 


6.5 


7.6 


(B) 


5.6 


6.9 


8.6 



TABLE I; The values of the ratio {2Ag/Tp} for all the his- 
tograms shown in Figs. fTHT b') and lllf bl. The first column 
corresponds to the mean-field results, while the second and 
third ones represent their shift under the influence of weak 
and strong fluctuations, respectively. The middle row is re- 
lated to model (A), and correspondingly to Fig. ITOT b). The 
bottom one contains the average ratios for model (B) and, 
consequently, for the histograms in Fig. Illf b') . 

fluctuations on the 2Ag/Tp distributions for models (A) 
and (B), respectively. The bottom histograms corre- 
spond to the mean-field samples, while the middle ones 
are related to the weaker fluctuations with 7 = 5 and the 
upper distributions are calculated for the more stronger 
fluctuations modelling by 7 = 10. It is clearly seen for 
the both models that the increase of the phase variation 
monotonously shifts the distribution to more higher val- 
ues without considerable changing of its shape. The av- 
erage values of the ratio 2Ag/Tp corresponding to all the 
histograms are represented in the Table HI It is worth to 
note that the average ratios corresponding to the stronger 
fluctuations with 7 = 10 (last column in Tab. |TJ are in 
excellent agreement with the appropriate ratios for the 
single perturbations. This points to the fact that these 
fluctuations suppress the proximity effect between neigh- 
bor enhanced pairing regions quite efficiently, so that the 
limit of single perturbation is practically reached for an 
individual perturbation region. 

The results demonstrated in Figs. [TOl E] and Tab. |T] 
are calculated for a given random realization of the pa- 
rameters {Oj}. Having studied a number of random re- 
alizations of these parameters we found that the absolute 
error in determining the average ratio (2Ag/Tp) is ~ 0.2. 

IV. ANTICORRELATION BETWEEN 
LOW-TEMPERATURE GAP AND 
HIGH-TEMPERATURE ZERO-BIAS 
CONDUCTANCE 

As it was shown above, one of the characteristic fea- 
tures of the model with the atomic-scale inhomogeneity of 
the superconducting OP is the considerable enhancement 
of the ratio 2Ag/Tp as compared to the homogeneous 
case. In this section we would like to discuss another 
characteristic manifestation of the OP inhomogeneity in 
the spectra measured by STM. It was recently reported'^^ 
that the value of the low-temperature position-dependent 
gap obtained from the STM spectra strongly anticorre- 
lated with the value of zero-bias differential conductance 
at a temperature, where practically the entire sample is 
ungapped. One of possible explanations for this experi- 
mental observation can be naturally given in the frame- 



work of the inhomogeneous OP model. 

Fig. [12] demonstrates the low-temperature gap-map 
(panel (a)) in comparison with the maps of high- 
temperature zero-bias thermally smeared LDOS (pan- 
els (b)-(d)). All the maps correspond to the considered 
above model (A) in the mean-field approximation. High- 
temperature zero-bias thermally smeared LDOS is cal- 
culated according to the formula ([4]) at V=0 for three 
different temperatures: the results for T — 0.15, where 
80% of the sample are still gapped, are shown in panel 
(b), while panel (c) is related to T = 0.19, where only 2% 
of the sample are gapped and panel (d) corresponds to 
T ~ 0.21, where the entire sample is already ungapped. 
The anticorrelation between the maps in panel (a) from 
one hand and the maps in panels (b)-(d) is clearly seen. 
Fig. fTST a) further quantifies this anticorrelation. High- 
temperature zero-bias thermally smeared LDOS, mea- 
sured along the vertical axis, versus low-temperature gap 
(horizontal axis) is plotted in this figure for all the sites 
of our sample. 

As it was described above in detail, the atomic-scale 
OP inhomogeneity leads to the fact, that the difference 
between the temperature of vanishing the gap Tp{i) and 
the temperature of disappearing the pair correlations, 
which corresponds to the critical temperature in the 
mean-field treatment, is considerably higher for the in- 
homogeneous situation as compared to the homogeneous 
case. The temperature evolution of thermally smeared 
LDOS from T = 0.15 to T = 0.21 is demonstrated in 
Figs. fTST bl-fd) for three typical space locations. Only 
very small part of the sample exhibits the behavior shown 
in panel (b), while approximately 2/3 of the sample can 
be described by panel (c) and the behavior of the con- 
ductance for the last 1/3 of the sites corresponds to panel 
(d). 

That is, although practically the entire sample is un- 
gapped at r = 0.19 for our model (A) and the conduc- 
tance curves change only slightly or practically do not 
change under further increasing in the temperature, the 
value of the superconducting OP is not small for the most 
part of the sample at this temperature (see Fig. El^a)). 
Even the complete vanishing of the gap at T = 0.21 does 
not mean that there are no pair correlations in the sys- 
tem. The presence and the strength of the pair corre- 
lations is reflected in the suppression of the low-energy 
LDOS. It is clearly seen in Figs. \W[a.)-{d) and \W[a.) 
that the anticorrelation weakens with temperature. It 
is rather weak at T = 0.21, where the pair correlations 
are small, and should disappear in the framework of our 
model when the pair correlations entirely vanish. We be- 
lieve that this scenario should also work if one takes into 
account the thermal phase fluctuations because the same 
physical picture is valid for the single perturbation as 
well. However, a regular theory is needed for a detailed 
consideration of this problem. 
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FIG. 12: (a) The low-temperature (T = 0.05) gap map for 
model (A), (b)-(d) The zero-bias space distributions of the 
thermally-smeared LDOS for the same model: (b) T — 0.15, 
(c) T = 0.19 and (d) T = 0.21. 

V. CONCLUSIONS 

We have studied the influence of the atomic-scale in- 
homogeneities of the superconducting OP on the conduc- 
tance spectra measured by STM. First of all, it is found 
that the ratio of the local low-temperature gap in differ- 
ential conductance spectra to the local temperature of 
vanishing the gap 2Ag/Tp can take large enough values 
as compared to the homogeneous OP model. While in 
the framework of the mean-field approximation the ratio 
does not strongly differ from the homogeneous one, the 
thermal phase fluctuations considerably enhance it. At 



least in the framework of a very simplified model we ob- 
tained that the ratio 2Ag/Tp can reach the values ~ 7—8, 
which are comparable to the experimental ones. It is also 
demonstrated that the additional weak potential scatter- 
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FIG. 13: (a) High-temperature zero-bias thermally-smeared 
LDOS, measured along the vertical axis, versus low- 
temperature gap (horizontal axis) for all the sites of the sam- 
ple at three different temperatures: T = 0.15, T — 0.19 and 
T = 0.21. (b)-(d) The temperature evolution of thermally 
smeared LDOS for three typical spatial locations. T = 0.15 
(solid black lines), T = 0.19 (dashed lines) and T = 0.21 
(dotted lines). All the panels are related to model (A). 



ers and hopping matrix element disorder do not influence 
qualitatively the above results. Second, the atomic-scale 
OP inhomogeneity results in the anticorrelation between 
the low-temperature gap and the high-temperature zero- 
bias conductance, bearing a resemblance to the recent 
resuhs obtained by STM. 
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